How generic are null spacetime singularities? 
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The spacetime singularities inside realistic black holes are sometimes thought to be spacelike 
and strong, since there is a generic class of solutions (BKL) to Einstein's equations with these 
properties. We show that null, weak singularities are also generic, in the following sense: there is a 
class of vacuum solutions containing null, weak singularities, depending on 8 arbitrary (up to some 
inequalities) analytic initial functions of 3 spatial coordinates. Since 8 arbitrary functions are needed 
(in the gauge used here) to span the generic solution, this class can be regarded as generic. 
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One of the most fascinating outcomes of general rela- 
tivity is the observation that the most fundamental con- 
cept in physics — the fabric of space and time — may 
become singular in certain circumstances. A series of sin- 
gularity theorems [jjj imply that spacetime singularities 
are expected to develop inside black holes. The observa- 
tional evidence at present is that black holes do exist in 
the Universe. The formation of spacetime singularities 
in the real world is thus almost inevitable. However, the 
singularity theorems tell us almost nothing about the na- 
ture and location of these singularities. Despite a variety 
of investigations, there is today still no consensus on the 
structure of singularities inside realistic black holes. 

At issue are the following features of singularities: their 
location, causal character (spacelike, timelike or null), 
and, most importantly, their strength. We use here 
Tipler's terminology for weak and strong singularities. 
In typical situations, if the spacetime can be extended 
through the singular hypersurface so that the metric ten- 
sor is C° and non-degenerate, then the singularity is weak 
|Q . The strength of the singularity has far-reaching phys- 
ical consequences. A physical object which moves to- 
wards a strong curvature singularity will be completely 
torn apart by the diverging tidal force, which causes un- 
bounded tidal distortion. On the other hand, if the sin- 
gularity is weak, the total tidal distortion may be finite 
(and even arbitrarily small), so that physical observers 
may possibly not be destroyed by the singularity [||,[| . 

The main difficulty in determining the structure of 
black hole singularities is that the celebrated exact black 
hole solutions (the Kerr-Newman family) do not give a 
realistic description of the geometry inside the horizon, 
although they do describe well the region outside. This 
is because the well known no-hair property of black holes 
— that arbitrary initial perturbations are harmlessly ra- 
diated away and do not qualitatively change the space- 
time structure — only applies to the exterior geometry. 
The geometry inside the black hole (near the singularity 
and/or the Cauchy horizon) is unstable to initial small 
perturbations [^|^] , and consequently we must go beyond 
the classic exact solutions to understand realistic black 
hole interiors. To determine the structure of generic sin- 
gularities, it is necessary to take initial data correspond- 
ing to the classic black hole solutions, make generic small 
perturbations to the initial data, and evolve forward in 
time to determine the nature of the resulting singularity. 
For this purpose a linear evolution of the perturbations 
may be insufficient — the real question is what happens 
in full nonlinear general relativity. 

The simplest black-hole solution, the Schwarzschild so- 
lution, contains a central singularity which is spacelike 
and strong. For many years, this Schwarzschild singu- 
larity was regarded as the archetype for a spacetime sin- 
gularity. Although this particular type of singularity is 
known today to be unstable to deviations from spherical 
symmetry (and hence unrealistic) ||, another type of a 
strong spacelike singularity, the so-called BKL singular- 
ity 0, is believed to be generic (below we shall further 



explain and discuss the concept of genericity). Since the 
BKL singularity is so far the only known type of generic 
singularity, in the last two decades it has been widely 
believed that the final state of a realistic gravitational 
collapse must be the strong, spacelike, oscillatory, BKL 
singularity. 

Recently, there have been a variety of indications that 
a spacetime singularity of a completely different type ac- 
tually forms inside realistic (rotating) black holes. In 
particular, this singularity is null and weak, rather than 
spacelike and strong. The first evidence for this new pic- 
ture came from the mass-inflation model H|| — a toy- 
model in which the Kerr background is modeled by the 
spherically-symmetric Reissner- Nordstrom solution, and 
the gravitational perturbations are modeled in terms of 
two crossflowing null fluids. Later more realistic analyses 
replaced the null fluids by a spherically symmetric scalar 
field 0. More direct evidence came from a non-linear 
perturbation analysis of the inner structure of rotating 
black holes jlO| . Both the mass-inflation models and the 
nonlinear perturbation analysis of Kerr strongly suggest 
that a null, weak, scalar-curvature singularity develops 
at the inner horizon of the background geometry. (See 
also an earlier model by Hiscock 

Despite the above compelling evidence, there still is a 
debate concerning the nature of generic black-hole singu- 
larities. It is sometimes argued that the Einstein equa- 
tions, due to their non-linearity, do not allow generic so- 
lutions with null curvature singularities p"2| . According 
to this argument, the non-linearity, combined with the 
diverging curvature, immediately catalyzes the transfor- 
mation of the null curvature singularity into a strong 
spacelike one — presumably the BKL singularity p3| . 
A similar argument was also given, some time ago, by 
Chandrasekhar and Hartle ||. According to this point 
of view, the results of the non-linear perturbation anal- 
ysis of Kerr are to be interpreted as an artifact of the 
perturbative approach used |]T5|| (and the mass-inflation 
model is a toy- model, after all). This objection clearly 
marks the need for a more rigorous, non-perturbative, 
mathematical analysis, to show that a generic null weak 
singularity is consistent with Einstein's equations. 

Recently, Brady and Chambers showed that a null sin- 
gularity could be consistent with the constraint section 
of Einstein's equations formulated on null hypersurfaces 
p6| . However, their result does not completely resolve 
the above issue. The hypothesis raised in Ref. |l2|, ac- 
cording to which nonlinear effects will immediately trans- 
form the singular initial data into a spacelike singular- 
ity, is not necessarily inconsistent with the analysis of 
Rcf . [[l6| . It is possible that a spacelike singularity could 
form just at the intersection point of the two charac- 
teristic null hypersurfaces considered in Ref. |T(|. It is 
primarily the evolution equations which will determine 
whether singular initial data will evolve into a null sin- 
gularity or into a spacelike one. 

The purpose of this paper is to present a new mathe- 
matical analysis which addresses the above question. Our 
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analysis shows that (i) the vacuum Einstein equations 
(both the constraint and evolution equations) admit so- 
lutions with a null weak singularity, and (ii) the class of 
such singular solutions is so large that it depends on the 
maximum possible number of independent functional de- 
grees of freedom. We will call such classes of solutions 
functionally generic (see below). Therefore any attempt 
to argue, on local grounds, that a null weak singularity 
is necessarily inconsistent with the non-linearities of Ein- 
stein's equations, must be false. In the present letter we 
outline this analysis and present the main results; a full 
account of this work is given in Ref . . 

Let us first explain what we mean by "degrees of free- 
dom" and "functionally generic" . Suppose that ip is some 
field on a 3+1 dimensional spacetime, which may be a 
multi-component field. Suppose that initial data for ip 
are specified on some spacelike hypersurface S. We shall 
say that ip has k "degrees of freedom" if k is the number 
of initial functions (i.e. functions of the 3 spacelike co- 
ordinates parameterizing S) which need to be specified 
on S in order to uniquely determine inside D + (S) the 
solution to the field equations satisfied by ip fig ). The 
number k depends on the type of field, and also possibly 
on the gauge condition used if there is gauge freedom. 
For example, for a scalar field k = 2, because one needs 
to specify both ip and ip on S. For the gravitational field, 
it is well known that there are 2x2 = 4 inherent degrees 
of freedom. The actual number k, however, is 4 plus the 
number of unfixed gauge degrees of freedom, which de- 
pends on the specific gauge conditions used. In the gauge 
we use, we find that k = 8 (see below). 

We shall say that a class of solutions to the field equa- 
tions is functionally generic, if this class depends on k 
arbitrary functions of three independent variables [fl9f . 
This concept of genericity is basically the same as that 
used by BKL [0. The motivation behind this defini- 
tion is obvious: Suppose that a given particular solution 
admits some specific feature (e.g. a singularity of some 
type). Obviously, in order for this feature to be stable to 
small (but generic) perturbations in the initial data, it is 
necessary that the class of solutions satisfying this fea- 
ture should depend on k arbitrary functions. Functional 
genericity is thus a necessary condition for stability, and 
is also necessary in order that there be an open set in 
the space of solutions with the desired feature, in any 
reasonable topology on the space of solutions [^(J . 

As we mentioned above, our result is a mathematical 
demonstration of the existence of a functionally generic 
null weak singularity. More specifically, we prove that 
there exists a class of solutions (M, g) to the vacuum 
Einstein equations, which all admit a weak curvature 
singularity on a null hypersurface, and which depend on 
k = 8 (see below) arbitrary analytic functions of three 
independent variables. (In Ref. |17|] we shall give a more 
precise formulation of this statement.) The singularities 
may also be characterized by the fact that the manifold 
may be extended through the null surface to an analytic 
manifold (M',g') where the metric g' is analytic every- 



where except on the null surface where it is only C°. 
Our construction is local in the sense that the manifolds 
we construct are extendible (in directions away from the 
null singularity); roughly speaking they can be thought 
of as open regions in a more complete spacetime, part of 
whose boundary consists of the singular null hypersur- 
face. We do not prove that null weak singularities arise 
in the maximal Cauchy evolution of any asymptotically 
fiat, smooth initial data set. The spacetimes we construct 
are of the form Z? + (E), where S is an open region in an 
analytic initial data set. The curvature singularity is al- 
ready present on the boundary of E in the initial data, in 
the sense that curvature invariants blow up along incom- 
plete geodesies. We emphasize that we do not view E as 
a physically-acceptable initial hypersurface; rather, the 
initial hypersurface E is merely a mathematical tool that 
we use to construct and parameterize the desired class of 
vacuum solutions. 

We shall first demonstrate the main idea behind our 
mathematical construction by applying it to a simpler 
problem — a scalar field. Consider, as an example, a 
real scalar field <p m flat spacetime, satisfying the non- 
linear field equation 

K = vtt) (i) 

where V(<f>) is some non-linear analytic function. (We 
add this non-linear piece in order to obtain a closer anal- 
ogy with the non-linear gravitational case.) In order to 
show that this field admits a functionally generic null 
singularity, we proceed as follows: Let x,y,u,v be the 
standard, double-null, Minkowski coordinates (i.e. such 
that ds 2 = —Adudv + dx 2 + dy 2 ). Equation (|l|) reads 

0,™ = ~ V(<l>) (2) 

where here and below the indices a, 6, ... run over the 
coordinates x and y. We now define 

W EE V 1/n (3) 

for some odd integer n > 3. We also define 

t EE w + u , z ee w — u. (4) 

Re-expressing the field equation (||) in terms of t and z, 
we obtain 

<P,tt = <P,zz + n[(z + t)j2] n - X - V (</>)] . (5) 

Let M° denote some neighborhood of the origin (x — 
y = z = t = 0) with compact closure, and let S + be the 
intersection of the hypersurface t=0 with M° (see Fig. 
1). Let fi(x,y,z) and f2(x,y,z) be two analytic func- 
tions of their arguments, defined on S + . For any such 
pair of functions, there exists a neig hborhood M+ C M° 
of S + , and a unique analytic solution cf>(x, y, z, t) to the 
field equation (J|) in M + , such that on S + , <p> = f\ 
and <f> t — f%- This follows directly from the Cauchy- 
Kowalewski theorem 12 ill , in view of the form of Eq. (g). 
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Let us denote the intersection of M + with the null hy- 
persurface v — by N + . Recall that N + includes a 
neighborhood of the origin in the hypersurface v = 0. 

Returning now to the original independent variables 
(u, v), we find that (j)(x,y,u,v) is continuous through- 
out M + . We now focus attention on the section v < 0, 
t > of M + , which we denote by M. Since the 
transformation from [z, t) to (u,v) is analytic as long 
as v ^ 0, we find that (f>(x,y,u,v) is analytic through- 
out M. However, (j> will generally fail to be smooth at 
v = 0: <f> v = {l/n)v 1 l n ~ 1 4>, w will diverge at v = as 
long as jto ^ there. We assume that at the origin, 
dfi/dz ^ if 2- This ensures that at least in some neigh- 
borhood of the origin, both cf> tW and cj> tU are nonzero. 
Let N be the intersection of that neighborhood with the 
section t > of N + . We find that (p >v diverges on N. 
Moreover, the invariant 4>. a <t>' a diverges on N too [it is 
dominated by (l/n)v~ 1+1 / n <j).u4 l .w ]• N is thus a singular 
null hypersurface. 
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FIG. 1. Spacetime diagram in z — t coordinates, illustrat- 
ing the mathematical construction used. Our final spacetime 
(M, g) consists of the shaded region. 



We conclude that there exists a class of solutions to 
Eq. (Q), which depends on two analytic functions of 
(x, y, z) (/i and f%) that can be chosen arbitrarily (apart 
from the above inequality), and which contains a singu- 
larity on a null hypersurface. In other words, the scalar 
field admits a functionally generic null singularity. (Note 
that 4> has a well- defined limit on the singular hypersur- 
face; this is the scalar-field analog of the notion of weak 
singularity.) 

We turn now to generalize this construction to the 
gravitational field. As before, our coordinates are de- 
noted (x,y, u, v). We adopt the gauge 

9ux 9uy Quu 9vv 0; (6) 

which in turn implies that g vv = 0. This ensures that the 
coordinate v is null (that is, the hypersurfaces v = const 
are null). There are six non-trivial metric functions, 
which we denote by gi (i = 1, •••,6), where here and 
below the indices i, j, . . . represent the six pairs of space- 
time coordinates (xx, xy, yy, vx, vy, uv). 



In this gauge, the number k of arbitrary functions in a 
general solution is k — 8. This can be seen as follows. De- 
fine the new variables T = v + u , Z = v — u. Then to de- 
termine a solution of the evolution equations, twelve ini- 
tial functions need to be specified on the spacelike hyper- 
surface T — const, namely gi(x,y,Z) and gi t r(x,y, Z), 
1 < i < 6. However, these 12 functions must satisfy 4 
constraint equations, as is always the case in general rel- 
ativity, so that the number of independently specifiable 
functions is k = 8. This conclusion can also be reached 
by adding the conventional number of intrinsic degrees 
of freedom of the vacuum gravitational field (2x2 = 4) 
to the number of unfixed gauge degrees of freedom in the 
gauge (§), which we show in Ref. |l7f to be 4. 

We shall now outline the generalization of the above 
scalar- field construction to the gravitational field. First, 
one writes the Einstein equations Rap = in the gauge 
(^). These equations can be naturally divided into six 
evolution equations and four constraint equations. At 
this stage we focus attention on the evolution equations, 
which can be taken to be Ri = 0. Next, we define w, 
t and z as before [Eqs. (H),(H)], and transform the field 
equations from the independent variables (u, v) to (z, t). 
[To avoid confusion, we emphasize that what we are do- 
ing here is not a coordinate transformation: it is just a 
change of independent variables in the differential equa- 
tions R a/ 3 = 0; thus, the unknowns in Eq. um below are 
still the six metric functions gi, which correspond to the 
coordinates (x,y,u,v).] By taking certain linear combi- 
nations of the equations Ri — 0, it is possible to rewrite 
the evolution equations in the schematic form 

9i,tt = fi(gj , 9j,t , 9j ,A , 9j,AB , 9j,At , z,t). (7) 

Here, the indices A, B run over the "spatial" variables 
x,y,z. If we impose certain inequalities on the ini- 
tial data [which ensure that in the region of interest 
det((/) = — 1], then the functions fa are analytic in all 
their arguments. [The gauge conditions (|^) are crucial in 
deriving Eq. (0).] 

We now consider the evolution of initial data under the 
system (^). As before, we take the initial hypersurface 
to be t — 0. Equation (^) requires twelve initial func- 
tions to be specified on this hypersurface: the six func- 
tions hi(x, y, z) = gi(x,y, z ; t = 0), and the six functions 
Pi(x,y,z) ee g i>t (x,y,z; t = 0). The form of Eq. (0) is 
suitable for an application of the Cauchy- Kowalewski 
theorem. Thus, defining S + , M + and M as before, and 
following the arguments above, we arrive at the following 
conclusion: For any choice of the above twelve analytic 
functions hi(x,y,z) and pi(x,y,z) on the section S + of 
t=0 (subject to certain inequalities), there exists an an- 
alytic solution gi(x,y,z,t) to Eq. (pi) in M. Again, re- 
turning from the variables (z,t) to the original indepen- 
dent variables (u,v), we find that the metric functions 
gi(x,y,u,v) are continuous throughout M + (and in par- 
ticular at v — 0) and, moreover, are analytic throughout 
M . However, at the hypersurface v = 0, gi <v typically 
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diverge like t> _1+1 /". As a consequence, the Riemann 
components R av bv generically diverge there jl7). More- 
over, it can be shown that the scalar K = R a p 1 sR a>3 ' yS 
also generically diverges at v — (like v~ 2+1 ^ n ). How- 
ever, it is easy to check directly that the singularity is 
weak. Thus, focusing attention on the physical region 
M , we find that the solutions constructed in that way 
are absolutely regular inside the region M, but develop 
a null, weak, scalar-curvature singularity on the portion 
v = of its boundary. 

The twelve initial functions hi(x,y,z) and Pi(x,y,z) 
are subject to four constraint equations. It is there- 
fore natural to expect that eight of these 12 initial func- 
tions can be chosen arbitrarily. This is not trivial to 
prove mathematically, however, especially because the 
constraint equations (expressed in the variables x, y, z) 
are somewhat pathological at z = 0. After some ef- 
fort, we found a mathematical construction which proves 
the above statement. More specifically, in our mathe- 
matical scheme one is free to choose the six hi(x,y,z), 
p xy (x,y, z), and one other function p(x,y,z). We can 
then show (using the Cauchy-Kowalewski theorem) the 
existence of a solution of the constraint equations (in a 
neighborhood of z = 0). The remaining initial functions 
Pi(x,y,z) are then determined from that solution. The 
above eight analytic functions can be chosen arbitrarily, 
up to some inequalities. 

To summarize, our mathematical construction shows 
the existence of a class of (local) solutions to the vac- 
uum Einstein equations, which contain a weak scalar- 
curvature singularity at the null hypersurface v — 0, and 
which depends on k = 8 analytic functions of (x,y,z). 
Our construction therefore demonstrates the existence of 
a functionally generic null, weak, scalar-curvature singu- 
larity. 

The main limitation of our construction is its restric- 
tion to analytic initial functions. We believe that this is 
merely a technical limitation of the mathematical theo- 
rems used in our proof, and the same physical situation 
(a null weak singularity) will evolve even if the initial 
functions on S + are not analytic (provided they are suf- 
ficiently smooth for v < 0). At any rate, it is worth- 
while to compare the mathematical status of our generic 
null weak singularity to that of the BKL singularity. To 
the best of the authors' knowledge, the existence of even 
a single inhomogeneous singular vacuum solution of the 
BKL type has not yet been proved mathematically - let 
alone the generality of this class of singular solutions. On 
the other hand we have demonstrated rigorously the ex- 
istence of a huge class of exact solutions containing null 
weak singularities. 

Our results have a simple intuitive interpretation: In 
linear hyperbolic systems, it is well known that weak dis- 
continuities of various types can freely propagate along 
characteristic lines. Our construction demonstrates that 
Einstein's equations, despite their nonlinearity, also be- 
have in this way (at least with respect to the type of weak 



discontinuity considered here). This is perhaps contrary 
to what was sometimes thought in the past, but is not 
really surprising, because, after all, Einstein's equations 
are quasi-linear. Thus, what we have shown is the lo- 
cal consistency of null weak singularities with Einstein's 
equations, despite the non-linearity of the latter. The 
important issue of the onset of the singularity from regu- 
lar, asymptotically flat, initial data (e.g. in gravitational 
collapse) still remains open; this issue is addressed (indi- 
rectly) by the nonlinear perturbation analysis of Ref. jl0| , 
but the onset still lacks a rigorous mathematical proof. 

Finally, it should be pointed out that the inner- 
horizon singularity that is suggested by perturbation 
analyses in a realistic rotating or charged black hole (see, 
e.g. Ref. JhJ), is qualitatively similar to the singularity 
constructed here, in that it is null and weak. There are 
some important differences, however. The main differ- 
ence is that the structure of the inner- horizon singularity 
is analogous to what would have been obtained from our 
construction if we had set w = l/ln|w|. Our method of 
proof does not generalize straightforwardly to this case, 
however, because v is no longer analytic in w at w = 
(though it is still C°°) p2fl . We hope to discuss the ana- 
lytic features of this more realistic null weak singularity 
elsewhere. 
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